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Abstract We study Coulomb drag between two parallel 
disordered mesoscopic ID-wires. By numerical ensemble av- 
eraging we calculate the statistical properties of the transcon- 
ductance G21 including its distribution. For wires with mu- 
tually uncorrelated disorder potentials we find that the mean 
value is finite, but with comparable fluctuations so that sign- 
reversal is possible. For identical disorder potentials the mean 
value and the fluctuations are enhanced compared to the case 
of uncorrelated disorder. 

1 Introduction 

Current flow in a conductor can through a Coulomb me- 
diated drag-force accelerate charge-carriers in a nearby 
conductor, thus inducing a drag-current. The effect is 
active whenever the distance between the two conduc- 
tors is of the same order as the distance between the 
charge-carriers - otherwise it is suppressed by screening. 
In the past years Coulomb drag in extended 2D-systems 
has been studied extensively [|] and very recently the 
study of fluctuations of the Coulomb drag was initiated 
by Narozhny and Aleiner || who found that the fluctu- 
ations will be pronounced for temperatures smaller than 
the Thouless energy. We study how drag between disor- 
dered mesoscopic ID-wires ]|,|| give rise to these new 
interesting phenomena such as a large fluctuations and 
sign reversal of the drag current. 

2 Formalism 

Consider two ID-wires of length L (shorter than the 
phase-breaking length) parallel to each other with a sep- 
aration d, see Fig. [p. Writing the Laplacian with the help 
of finite differences the Hamiltonian of the uncoupled 
wires is mapped onto a tight-binding model 
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where i = 1, 2 label the two wires and n, n' = 1, 2, 3, . . . N 
the lattice points. The conducting properties can be ob- 
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Fig. 1 Coulomb coupled ID-wires where • denote the lattice 
points of the wires and o denote those belonging to the ideal 
leads. 



taincd from the retarded Green functions of the iso- 
lated wires which can be written as N x N matrices: 
Qi = (ep — Hi — S % L — Z 1 ^) -1 , where S l p is the retarded 
self-energy describing coupling to the lead p — L, R. 

Using Kubo formalism we calculate the Coulomb drag 
to second order in the interaction U\2 between the meso- 
scopic ID-wires which we assume to be otherwise non- 
interacting ||||] . For kT <^ s F the dc transconductance 
G21 = dh/dVi becomes §,§ 

G21 = j (kTf jTt[U 12 M 1 Uu M a ] , (2a) 

where U12 is an N x N coupling matrix representing 
the interwire Coulomb interaction and Mj is an TV x N 
matrix 



(2b) 



Mi = Rc{A[ <g> [AiAAi] } , Ai = i [G, - g\] . 

Here, A nn , = ±5 n>n ,±x/{N - 1) and {X <g> Y} nn , = 
X nn 'Y nn >. The Landauer conductance Gu — dh/dVi of 
the individual wires can be expressed in a similar form Q| 



g u = ^iv[r| g i r* R g\] ,r; = i[s;- {s;} 



(3) 



3 Ensemble averaging 

The statistical properties of drag can be analyzed by 
generating an ensemble of different disorder configura- 
tions and using Eq. (0) to calculate the drag. For the 
disorder we use the Anderson model with diagonal dis- 
order H where the transport mean free path £ can be 
related to the disorder strength W by I = al2(4tej? — 
e F )/W 2 . We consider two cases: i) both wires being dis- 
ordered, but with U\ and U2 fully uncorrelated and ii) 
both wires being disordered and fully correlated ||, i.e. 
U\ = 1/2- For weak disorder a diagrammatic perturba- 
tion expansion for the fluctuations SG21 = G21 — (G21) 
gives ([«?2i] 2 ) 1/2 cx 1/kpt §| and it can also be ar- 

gued that ([£G 2 i] 2 ) c = 2 x <[<5G 2 i] 2 ) uc §. Both pre- 
dictions are valid to lowest order in l/kp£. 

4 Results 

We consider quarter-filled bands (ep = t) and wires with 
N = 100 lattice points so that k F L = (tt/3) x 100. 
The separation is kpd = 1 and for simplicity we as- 
sume an unscreened coupling of the form {EA.2} / = 
e 2 /(47re e r [(n-n') 2 a 2 +d 2 ] 1/2 ). 
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Fig. 2 Right panel: fluctuations (<5G|i {£)) 1/2 
system parameters). Left panels: histograms (based on 
k F l = (tt/3) x 3600. 



normalized by the ballistic result G2i(oo) as a function of k F £ (see text for 
10 4 random disorder configurations) for Gi\(€) in the case of 



We study the de-localized regime £ ^> L where we 
as expected |p!o[[riJ| find that the fluctuations 5Gu of the 
Landauer conductance Ga are vanishing and (Gjj) ~ 
2e 2 /h. However, for the transconductance G21 even weak 
disorder can have a large effect. The lower left panel 
of Fig. H shows a typical histogram of G2i(^)/G ! 2i(oo) 
(where G2i(oo) is the result in the ballistic regime, U% — 
U2 = 0) for £ — 36L. Depending on the disorder con- 
figuration G2i(^) can be either higher or lower than in 
the ballistic regime. The enhancement occurring for cer- 
tain disorder configurations can be understood physi- 
cally as follows. The lack of translational invariance al- 
lows forward scattering (transferred momentum q ~ 0), 
which normally has little effect, to cause transitions be- 
tween scattering states with opposite directions, thus 
contributing to the drag. The variance is of the same 
order as the mean value so that sign reversal for some 
disorder realizations is possible. The latter is represented 
by the negative tail in the histogram. 

For the same system parameters but now with iden- 
tical (correlated) disorder potentials we get a very differ- 
ent distribution as seen in the upper left panel of Fig. ||. 
As predicted |)| the mean value is enhanced compared 
to uncorrelated disorder and also the fluctuations are 
enhanced. In the right panel of Fig. |^ we show the de- 
pendence of the fluctuations on the mean free path fcpf 
which has the expected l/kp£ behavior. Comparing the 
two disorder situations we find numerical support for the 
predicted relative strength of y/2 H . 

5 Conclusion 

We have numerically studied drag of disordered meso- 
scopic ID-wires in the de-localized regime £ ^> L. Our 



results illustrate how the statistics of the transconduc- 
tance depend strongly on disorder and we find that even 
weak disorder can give rise to fluctuations of the same 
order of magnitude as the transconductance for the bal- 
listic case. This implies that the direction of drag de- 
pends on the disorder configuration and that for a given 
system the sign of the drag current will be arbitrary. Our 
results also confirm the for 2D extended systems recently 
predicted enhancement of the mean value for correlated 
disorder compared to uncorrelated disorder. In addition 
we have also found a corresponding enhancement of the 
fluctuations by a factor of \[2 compared to uncorrelated 
disorder. 

We acknowledge C. W. J. Beenakker and M. Brand- 
by ge for useful discussions. 
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